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ABSTRACT 

We determine the average metallicities of the eleraents of cold halo substructure (ECHOS) that we 
previously identified in the inner halo of the Milky Way within 17.5 kpc of the Sun. As a population, 
we find that stars kinematically associated with ECHOS are chemically distinct from the background 
kinematically smooth inner halo stellar population along the same Sloan Extension for Galactic Under- 
standing and Exploration (SEGUE) line of sight. ECHOS are systematically more iron-rich, but less 
a-enhanced than the kinematically-smooth component of the inner halo. ECHOS are also chemically 
distinct from other Milky Way components: more iron-poor than typical thick-disk stars and both 
more iron-poor and a-enhanced than typical thin-disk stars. In addition, the radial velocity disper- 
sion distribution of ECHOS extends beyond tr 20 km s~^. Globular clusters are unlikely ECHOS 
progenitors, as ECHOS have large velocity dispersions and are found in a region of the Galaxy in 
which iron-rich globular clusters are very rare. Likewise, the chemical composition of stars in ECHOS 
do not match predictions for stars formed in the Milky Way and subsequently scattered into the inner 
halo. Dwarf spheroidal (dSph) galaxies are possible ECHOS progenitors, and if ECHOS arc formed 
through the tidal disruption of one or more dSph galaxies, the typical ECHOS [Fe/H] ^ — 1.0 and 
radial velocity dispersion ct ~ 20 km s^^ implies a dSph with Mtot ^ 10^ Mq. Our observations con- 
firm the predictions of theoretical models of Milky Way halo formation that suggest that prominent 
substructures are likely to be metal-rich, and our result implies that the most likely metallicity for a 
recently accreted star currently in the inner halo is [Fe/H] ^ —1.0. 

Subject headings: Galaxy; abundances — Galaxy: formation — Galaxy: halo — Galaxy: kinematics 
and dynamics 



1. INTRODUCTION 

In ACDM cosmology, galaxies like the Milky Way are 
formed through the stochastic accretion of smaller dark 
matter halos and the continuous accretion o f gas (e.g., 
IPress fc Schechte^flOTl fWhite fc ReesI 119781 ). The sta- 
tistical properties of the stochastic accretion have been 
well characterized by large cosmological dark matter 
only n-body simulations (e.g., Springcl ct al. 2005). On 
the other hand, no model has yet self-consistently in- 
corporated the gas dynamics, star formation, and feed- 
back necessary to satisfactorily reproduce the formation 
of Milky Way analogs. In any case, an understand- 
ing of the statistical properties of the accretion histo- 
ries of Milky Way analogs does not predict the unique 
formation history of the Milky Way. Fortunately, the 
dynamical and chemical properties of the Milky Way's 
thin disk, thick disk, bulge, and halo provide strong 
constraints on the formation of our own Galaxy (e.g.. 
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iFreeman fc Bland-Hawthornll2002l : lHelmill2008| ). 

The kinematics of the stellar halo of the Milky Way 
are simpler to describe than the kinematics of other 
Galactic components. As a result, the stellar halo is the 
Galactic component in which it is easiest to unambigu- 
ously identify the remnants of the smaller halos that the 
Milky Way has accreted through its history. For that 
reason, the substructures recently discovered in the halo 
provide a direct measure of accretion in the region of the 
Galaxy that is far easier to model than the disk or the 
bulge. The existence of substructure in the halo of the 
Milky Way is n ow well established by star co unts (e.g., 
Totten fc IrwinI IT99 8: Tottcn et al. '2000; I vezic et alJ 
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iGUmore et al.l 120021: iNewberg et al l 
120021; iRockosi et'al.' '2002^ 'Maiewski et al.' '20031 
Yanny_eial. 2003; Rocha-Pinto ct al. 2004; Duffau ct al| 
2006s Belok urov et al.l 120061: iGriUmair fc Johnso: 
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Belokurov e t al. 2007; B ell et al 
Grillmaid 12009 : Watki ns et al 
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1996; Chiba & Yoslil 
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Keplev etal.1 120071; llyezic et al.1120081: iKlement et al l 



2008:: ISeabroke et al.l 120081: iKlement et al.l 120091: 
ISmith et al. l 120091: iStarkenburg et al.ll2009D. and chemi- 
cal abundances (e.g.. Ilvezic et al! I I2008I : lAn et"all [2009): 

these observed substructures are likely the remains of the 
stellar populations that formed as part of independent, 
bound structures that are now being disrupted and 
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dispersed in the halo. 

In t he fi rst paper in this series, iSchlaufman et afl 
(|2009[ )- [S09l hereafter-we identified elements of cold halo 
substructure (ECHOS) as overdensities in the spatial and 
radial velocity distribution of the inner halo's stellar pop- 
ulation. Except for those ECHOS identified along lines of 
sight known to host surface brightness substructure, none 
of our ECHOS are detected as overdensities of stars in 
surface brightness substructure. As a result, they are dis- 
tinct from surface brightness substructure and are likely 
the debris of more ancient accretion events. We used 
our detections to infer that as much as 10% of the inner 
halo by volume might have 30% of its stars in ECHOS, 
and we used that measurement together with similar 
measurements for surface brightness substructure to sug- 
gest that the Milky Way's accretion history has been 
roughly constant over the last few Gyr with no massive 
(Mtot ~ 10^° ^^o) single accretion events in that inter- 
val. Moreover, the distribution of ECHOS in Galactic co- 
ordinates is consistent with isotropy given the complete- 
ness limits of our search, possibly indicating an accretion 
origin for ECHOS (rather than an association with the 
Galactic disk). 

The chemical composition of the Milky Way's stellar 
popu lations also illuminate their origin (e.g.. iMcWilliamI 
Il997f) . In general, since iron is introduced into the inter- 
stellar medium in supernovae (SN) explosions, the [Fe/H] 
of a stellar population is correlated with the total inte- 
grated star formation that occurred in that population. 
Qualitatively then, low [Fe/H] indicates relatively little 
integrated star formation, while high [Fe/H] indicates rel- 
atively more integrated star formation. If the initial mass 
function (IMF) of stellar populations is not a strong func- 
tion of env ironment or metallicity at [Fe/H] > —3.0 (e.g., 
iBastianet al. 2010), then a stellar population's abun- 
dances of the a-elements (O, Mg, Si, Ca, and Ti) relative 
to iron are related to the duration of the star formation 
that produced that population. Core-collapse (i.e.. Type 
II and Type Ibc) SN result from the explosions of stars 
initially more massive than about 8 Mq and produce 
large amounts of the a-elements relative to iron. Since 
stars that massive spend very little time on the main se- 
quence, the enrichment of the interstellar medium (ISM) 
by the ejecta from core-collapse SN occurs very quickly, 
within a few Myr. In contrast, thermonuclear (i.e.. Type 
la) SN result from the explosion of white dwarfs as a 
result of runaway nuclear burning and produce large 
amounts of iron relative to the a-elements. The exact 
progenitor and therefore the characteristic timescale on 
which thermo nuclear SN start to enrich the ISM is still 
debated (e.g.. iScannapieco fc Bildstenl 120051 ): neverthe- 
less, the timescale is longer than the comparable time 
to ISM enrichment through core-collapse SN. For those 
reasons, a short burst of star formation will leave behind 
a stellar population with [a/Fe] > 0, while a more ex- 
tended star formation history will produce a population 
with [a/Fe] ^ (neglecting any differences in the IMF). 
Accordingly, a lower star formation rate will allow ther- 
monuclear SN to reduce [a/Fe] at relatively low [Fe/H], 
whereas a high star formation rate will produce many 
generations of stars before thermonuclear SN begin to 
reduce [a/Fe], pushing the "knee" in the [Fe/H]-[a/Fe] 
plane to higher [Fe/H]. 

Stars in the inner halo typically have [Fe/H] ^ —1.6 



and are enhanced in the g -el ements relative to iron 
(e.g.. IRvan fc Norrid Il991allbl: iMcWiUiam et al.l [l995l: 
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lAllende Prieto et al.H2006f) . [Robertson et al.l (|2005D ex- 
plained the chemistry of most of the stellar mass in 
the inner halo in the context of the ACDM model of 
galaxy formation with the accretion of a massive M^ot 
5 X 10^" Mq halo ~ 10 Gyr in the past. The high-mass 
and short-timescale for star formation in such a massive 
progenitor of the inner halo are consistent with the ob- 
served chemistry. This scenario is in contrast to the com- 
position of the classical dwarf spheroidal galaxies (dSph) , 
which at the averag e [Fe/H] of th e inner halo have [a/Fe] 
closer to solar (e . g..lMateo l ll998f) . At the same time, the 
[Robertson et al.l (|2005f) niodel did not directly address 
the origin of the substructure now known to exist in the 
inner halo, and those substructures are likely related to 
more recent accretion events. 

In the same way, the chemical properties of substruc- 
ture reveal something about its origin. Since ECHOS 
are likely the remains of recent accretion events (t < 5 
Gyr), their chemical composition constrains the prop- 
erties of a more recently-disrupted class of contribu- 
tors to the stellar population of the inner halo than 
the few massive progenitors that contributed the bulk 
of the inner halo's stellar population. If ECHOS are 
the remains of disrupted dSph galaxies, then the dSph 
luminosity-metallicity relation can be used to infer the 
luminosity of a progenitor fr om its mean iron metallicity 
(e.g iHelm i et al. 2006bi; iKirbv et al.l l2008bD . Indeed, 
the IKirbv et al.. (2008bh luminosity-metallicity relation 
for dSph galaxies is now well calibrated over four decades 
in lumino sity using. C ombined with stellar mass to light 
ratios, the lKirbv et al.l (l2008b) luminosity-metallicity re- 
lation can be turned into a stellar mass-metallicity rela- 
tion. As a result, the average iron metallicity of a dSph 
stellar population can be considered a proxy for its stel- 
lar mass. Likewise, the [a/Fe] of substructure indicates 
the duration of star formation in its progenitor. Typical 
halo stars have [a/Fe] > 0.3 at [Fe/H] ^ —1.6 indicating 
formation in a truncated episode of star formation. In 
contrast, stars in classical dSph typical have [a/Fe] ~ 
at [Fe/H] ~ —1.6, indicating formation in an episode of 
prolonged star formation. 

We measu re th e chemical properties of the ECHOS 
discovered in IS09I and use their combined kinematic and 
dynamical properties to better understand their origin 
in the context of the ACDM model of Galaxy forma- 
tion. This paper is organized as follows: in Section 2 we 
describe the data we use in this analysis. In Section 3 
we describe how we derive average [Fe/H] and [a/Fe] for 
each ECHOS through the analysis of coadded metal-poor 
main sequence turnoff star spectra creating using spectra 
from individual metal-poor main sequence turnoff stars 
kinematically associated with each ECHOS. In Section 4 
we discuss the implications of our findings for the forma- 
tion of the Milky Way. We summarize our conclusions in 
Section 5. 

2. DATA 

The Sloan Extension for Galactic Understanding and 
Exploration (SEGUE) survey observed approximately 
240,000 Milky Way stars with apparent magnitudes in 
the range 14 < g < 20.3 with the fiber-fed Sloan Digital 
Sky Survey (SDSS) spectrograph at moderate resolution. 
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Spectroscopic targets were selected from the combined 
11,663 deg^ ugriz photometric footprint of the SDSS 
and SEGUE. The SDSS telescope and spectrograph 
obtain R « 1800 spectra between 3900 A and 9000 
A with high spectrophotometric accuracy. The SEGUE 
instrumentation, data processing pipelines, survey 
strategy, along with radial velocity and atmospheric 
param ete r accuraci es are descri bed in [ Yannv et aT 
20091 ) ■ iLee et al.l (2008a,b), lAUende Prieto et al 



^08!), and the SDSS-II DR7 paper (jAbazaiian et al 
20091). The S PSS su rvey is described in detail 
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In IS09l we examined a sample of 10,739 metal-poor 
main sequence turnoff stars collected from 137 seven deg^ 
lines of sight. These metal-poor main sequence turnoff 
(MPMSTO) stars have both the g — r color and signif- 
icant UV excess expected for the main sequence turnoff 
of a metal-poor population (for a detaile d de scription 
of the MPMSTO sample see Section 2 of [SOl). Given 
the magnitude limits of the SEGUE survey, the MPM- 
STO sample was selected because MPMSTO stars are 
the highest-density tracer of the inner halo. In IS09I . we 
defined the inner halo as the volume more than 10 kpc 
from the Galactic center, within 17.5 kpc of the Sun, and 
more than 4 kpc from the Galactic plane. Though we 
found no reason to reject a kinematically smooth model 
for the inner halo on average, we discovered many radial 
velocity overdensities that we termed ECHOS. We iden- 
tified ECHOS along individual lines of sight along which 
there was a very significant deviation from a kinemati- 
cally smooth model of the inner halo. We gave these sub- 
structure the name ECHOS to differentiate them from 
surface brightness substructure like tidal streams, as the 
only ECHOS we could relate to surface brightness sub- 
structures were those ECHOS discovered along lines of 
sight targeted at known surface brightness substructures. 
ECHOS are also distinct in that they are likely related 
to more ancient accretion events than su rface brightness 
subst ructure in the same volume (e.g., iJohnston et al.l 
I2008f) . As the same time, ECHOS are likely related to 
more recent accretion events than those substructures 
that will be discovered in the same volume using the 6D 
kinematic data that will become available from Gaia and 
the LSST ('e.g.. IHelmi etlll lMifial: i McMiUan fc BinnevI 
[200l iGomez fc Helmil l2010t ). The ECHOS we discov- 
ered were preferentially bunched at the faintest appar- 
ent magnitudes of the SEGUE spectroscopic sample; in- 
deed the 25%, 50%, and 75% percentile in r-magnitude 
and spectral signal to noise per 1 A pixel (abbreviated 
S/N hereafter) were (19.0,19.5,19.8) and (7.6,10.4,15.4) 
respectively. As a result, the spectra of individual MPM- 
STO stars were generally at too low a S/N to precisely 
measure abundances. 

3. ANALYSIS 

We coadd all stellar spectra kinematically associated 
with an ECHOS in a narrow range of effective tempera- 
ture and analyze the resultant coadded spectrum to de- 
termine the average abundance of the ECHOS. At the 
same time, we use an equivalent coaddition process to 



determine the average metallicity of the MPMSTO stars 
in the kinematically smooth inner halo population along 
the same line of sight. 

We use these measurements to compare the mean 
chemical abundance of the ECHOS and smooth halo 
population. Both calculations are subject to the same 
systematics, so any apparent chemical offset can only re- 
sult from an underlying chemical difference or from ran- 
dom effects. To quantify the effects of randomness, we 
compute the precision and accuracy of our technique by 
analyzing two classes of objects with known composi- 
tion. The first class of objects are individual MPMSTO 
stars that have been observed both at high S/N during 
the SEGUE survey and at high resolution by larger tele- 
scopes. The high S/N SEGUE spectra can be degraded 
to arbitrarily low S/N by a detailed noise model so that 
we can test the coaddition process with data represen- 
tative of the low S/N spectra available for MPMSTO 
stars in ECHOS. The second class of objects are MPM- 
STO stars associated with the globular clusters M 13 and 
M 15. Using those two test cases, we quantify the mean 
square error (MSE = bias^ -I- variance) of our metallic- 
ity estimates. We describe these steps in detail in the 
following subsections. 

3.1. The SEGUE Stellar Parameter Pipeline 

The SEGUE Stellar Parameter Pipehne (SSPP - 
ILee et al. 2008a, b; AUcnde Prieto ct al. 2008) uses Sloan 
spectroscopy and ugriz photometry to infer the stellar 
atmosphere parameters (Toff, logg, [Fe/H], and [a/Fe]) 
of stars observed in the course of SDSS and SEGUE. 
The SSPP implements a multimethod algorithm in which 
many different techniques are used to compute the stellar 
parameters. The SSPP then averages the result of each 
method known to be valid in a given color and S/N range 
to determine the final Toff, logg, [Fe/H], and [a/Fe] re- 
ported for all stars observed in the SDSS and SEGUE 
s urveys. 

ILee et all (|2008ai rbl) determined the accuracy and pre- 
cision of the SSPP in three ways. First, they compared 
the atmospheric parameters determined by the SSPP 
from high S/N SEGUE spectra with the atmospheric pa- 
rameters determined from high-resolution spectroscopy 
from HIRES and ESI on the Keck Telescopes, HRS on 
the Hobby-Eberly Telescope (HET), and HDS on the 
Subaru Telescope. Their results suggest that including 
both systematic and random error, the SSPP has a one- 
sigma precision of 141 K in T^s, 0.23 dex in logg, and 

0. 23 dex in [Fe/H]. The stars bright enough to be ob- 
served at high resolution all had S/N > 50 SEGUE spec- 
tra, so direct comparison with high-resolution spectra 
can only be made for high S/N spectra. For that rea- 
son, they degraded these high S/N SEGUE spectra with 
stellar parameters well characterized by high-resolution 
spectroscopy with a detailed noise model to create many 
thousand spectra at many values of S /N between 55 and 

1. They then ran the SSPP on the noise-degraded spec- 
tra to determine the accuracy and precision of the SSPP 
as a function of S/N and reported the result in in Table 
6 of ILee et all (|2008aD . They found a [Fe/H] precision of 
(0.5,0.2,0.1) dex at S/N of (5,10,15) andaJa/F^ preci- 
sion of 0.1 dex at S/N gre ater than 20 ()Lee et al.ll2011l) . 
Finally, ILee et al.l ()2008l3[ ) applied the SSPP to stars as- 
sociated with open and globular clusters, where they 
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found that the SSPP achieves a precision in [Fe/H] of 0.13 
dex over the range —0.3 < g ~r < 1.3, 2.0 < logy < 5.0, 
and -2.3 < [Fe/H] < 0.0. 

3.2. Coaddition Algorithm 

We include in the coadded spectra only those spectra 
that correspond to MPMSTO stars with radial velocities 
that indicate membership in the population of interest. 
In addition, we include in the coadded spectra only those 
spectra that correspond to MPMSTO stars in a finite 
range of effective temperature, as the strength of spec- 
tral lines is affected by temperature as well as [Fe/H]. 
Consequently, the Tob range has to be small enough to 
ensure that the coadd gives an accurate [Fe/H] that is 
representative of the population. In order to estimate 
the effective temperature of a MPMSTO star from its 
g — r color, we fit a linear model between g — r color and 
Teff in the range 0.15 < g — r < 0.5 for all stars in SDSS 
DR7 with SEGUE spectroscopy at S/N > 40 and a SSPP 
[Fe/H] < -1.0. We find that Toff « -3800(5 " + 7300 
with about 500 K of scatter at constant g — r color. We 
use g — r color because the photometric accuracy of the 
SDSS does not vary over the apparent magnitude range 
of our MPMSTO sample, while the spectral S/N and 
therefore T^g precision inferred from the spectra varies 
substantially. Moreover, there are no reliable Teff esti- 
mates for the half our MPMSTO sample with S/N < 10. 
For those reasons, we use the simple relationship between 
g — r color and Tcs at high S/N to ensure that we reli- 
ably include in each coadd spectrum only MPMSTO star 
spectra in a narrow range of Toff even at low S/N. Like- 
wise, the photometric selection and spectroscopic con- 
firmation of MPMSTO stars indicates that the stars in 
our sample have similar surface gravities. As a result, 
we ensure that each spectrum that enters into a coadd 
correspond to a star in a narrow range of Tos and logg. 
Therefore, the only unconstrained stellar parameter is 
metallicity, and a SSPP analysis of the coadded MPM- 
STO spectrum will produce an unbiased estimate of the 
mean metallicity of the MPMSTO population. 

We shift each MPMSTO spectrum eligible for inclu- 
sion in a coadd to a heliocentric radial velocity Vr — 
km s^^. We then use natural cubic spline interpolation 
to interpolate both the spectrum and its inverse variance 
on to a common grid in wavelength. Next, we numeri- 
cally integrate the area under the curve defined by the 
spectrum and normalize both the spectrum (by divid- 
ing by the normalization factor) and the inverse variance 
(by multiplying by the normalization factor squared) to 
ensure that each spectrum that is to be included in the 
coadd has the same scale. For each population of in- 
terest, we then create an ensemble of realizations of the 
coadded spectrum by bootstrap resampling from the set 
of radial velocity zeroed, interpolated, and normalized 
spectra that belong to that population. Each spectrum 
contributes to each wavelength bin in proportion to its in- 
verse variance in that bin relative to the other spectra se- 
lected for coaddition. One danger to this approach is the 
possibility that the resultant coadded spectrum does not 
correspond to the spectrum of any physical star. This is 
unlikely in our analysis though, as we obtain good agree- 
ment between globular cluster metallicities produced by 
coaddition and their known metallicities from high reso- 
lution spectroscopy. We describe the coaddition process 



in detail in Appendix A. 

The SSPP also uses Sloan ugriz photometry in its pa- 
rameter estimation routines, so we determine an equiva- 
lent ugriz photometric measurement for our coadd spec- 
tra by computing a weighted average of the ugriz photo- 
metric measurement of the individual stars in each boot- 
strap coadd, using the mean S/N between 3950 A and 
6000 A as the weight. We then run the SSPP on each of 
the bootstrap coadded spectra and average ugriz pho- 
tometry to determine the mean of the SSPP [Fe/H] and 
[a/Fe] estimates for that particular MPMSTO popula- 
tion. 

3.3. Accuracy and Precision in the Ideal Case 

To determine the precision and accuracy of our coad- 
dition algorithm as a function of S/N and population 
metallicity, we analyzed spectra created by coadding in- 
dividual MPMSTO star spectra that had been observed 
to very high S /N during the SEGUE survey but that had 
been subsequently degraded with a detailed noise model. 
Though this is the ideal case of uniform effective tem- 
perature, surface gravity, and metallicity, the result will 
reveal the amount of bias and variance in our measure- 
ments that can be attributed to noise and population 
metallicity. 

Section 6 of iLee et al.l ()2008aD describes the algorithm 
used to create the noise degraded spectra we use in this 
analysis and summarize here. Each SEGUE plate ob- 
tains spectra for stars that span a range of magnitudes 
from g fa 15 to g 20, but is exposed to a fiducial S/N 
for the faint targets. As such, the noise properties of the 
SEGUE spectra vary with the magnitude and color of 
the targets, with increasing fractional contribution from 
the sky for faint targets. Because the observing crite- 
ria for the survey were homogeneous, it is possible to 
parametrize those variations and create a realistic model 
noise spectrum that can be used to create low S/N real- 
izations of high S/N spectra. Randomly chosen residual 
spectra from the sky fibers are then added to complete 
each realization. This noise model was used in ILee et al.l 
(jlOMb.) to test the accuracy of the SSPP, and we use the 
same model here to test the accuracy of the SSPP param- 
eters of our coadds. We use a sample of 640 noise-added 
reahzations of eight high S/N SEGUE MPMSTO stars 
with [Fe/H] values from -2.41 < [Fe/H] < -0.23, which 
spans the metallicity range of our sample. We have 54 
realizations at S/N values from 7.5 to 55. 

We use this large sample of noise-degraded spectra for 
eight MPMSTO stars to determine the precision and ac- 
curacy of our coaddition algorithm as a function of metal- 
licity and S/N. The analysis of noise-added spectra al- 
lows us to examine the performance of the SSPP and 
our coaddition algorithm over a range of metallicity and 
S/N that spans our entire sample. We plot the results 
of our analyses for [Fe/H] in Figured] and for [a/Fe] in 
Figure [H The MSE of our SSPP analysis of coadded 
spectra created from the coaddition of noise-degraded 
spectra of single MPMSTO stars ranges from 0.05 dex 
in both [Fe/H] and [a/Fe] for the most iron- rich stars to 
0.20 dex in both [Fe/H] and [a/Fe] for the most iron-poor 
stars. In both cases, the range in effective S/N results be- 
cause the coadded spectra are created by coadding from 
an ensemble of spectra with different S/N, and the range 
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of S/N apparent in Figure [T] and Figure [5] is larger than 
the equivalent range of S/N in the ECHOS. Reassuringly, 
there is no obvious trend in the precision and accuracy of 
the bootstrap coaddition process with effective tempera- 
ture. The outliers at low S/N in Figure[T]and Figure[2]are 
likely coadded spectra that include very low S/N spec- 
tra with large sky residuals, possibly taken from SEGUE 
plates with below-average sky subtraction. The existence 
of outliers due to this effect is another reason why the 
bootstrap resampling from many spectra is important 
to determine the error distribution that results from the 
coaddition process. In the next subsection, we will use 
this analysis of the noise-added data to determine the 
MSE of our coadditon analysis as function of metallicity 
and S/N. 

3.4. Accuracy and Precision in a Representative Case 

In reality, the spectra in each ECHOS coadd corre- 
spond to stars in a small but finite range of effective 
temperature, surface gravity, and composition. In order 
to determine the degree to which these spreads in prop- 
erties affects our ability measure the mean metallicity of 
the ECHOS and smooth halo populations by coadding 
spectra belonging to each, we analyzed spectra created 
by coadding individual MPMSTO spectra corresponding 
to MPMSTO stars that likely belong to the well studied 
Galactic globular clusters M 13 and M 15. The M 13 and 
M 15 data are representative of the range in Tcs and log g 
of a real MPMSTO population as selected by the SEGUE 
survey. It is exactly the fact that the globular cluster 
MPMSTO stars occupy a narrow range in [Fe/H] that 
makes the this globular cluster data so useful-we know 
the expected [Fe/H] very well and can test whether or 
not we converge to the known value when coadding from 
a range of Toff and logg. Unfortunately, [Fe/H] and Tog 
are degenerate, so we optimize the g — r range selected 
in order to produce coadded spectra with the maximum 
S/N and therefore the most accurate abundance estimate 
from the SSPP in the following way. Increasing the T^e 
range of spectra that are coadded increases the number 
of spectra included in each coadd, increases the total 
signal in each coadd, and ultimately increases the pre- 
cision of the abundance estimate. The trade off is that 
coadding spectra in too large a range in Tcs can decrease 
the accuracy of the abundance estimate. We find that 
coadding spectra that belong to stars in 250 K bins pro- 
duces the same accuracy as the coaddition of stars in 500 
K bins, as the scatter in the Tee—g — r relation at con- 
stant g — r color is about 500 K. Combined with the fact 
the coadd spectrum that results from coadding spectra in 
the 500 K bin always reaches higher S/N (and therefore 
higher precision) than the coadd spectrum that results 
from coadding spectra in the 250 K bin, we exclusively 
use the 500 K bin in our analysis of ECHOS. 

We select those stars that have equatorial coordinates 
that place t hem within t he tidal radius of each cluster (as 
reported in lHarrislll996f) and that have radial velocities 
consistent with cluster membership. Given the precision 
of SEGUE radial velocities at the S /N of the cluster spec- 
tra for cluster members that meet the SEGUE turnoff 
sample criteria, that corresponds to 15 km s~^ for M 13 
and 25 km s^^ for M 15. We select for coaddition only 
those spectra that correspond to MPMSTO stars with 
g — r colors that place them within the 500 K bin in effec- 



tive temperature that produces the highest S/N coadd. 
The median S/N of the MPMSTO spectra that remain 
after applying these cuts is 22.9 from 12 spectra for M 13 
and 9.0 from 11 spectra for M 15. The number of spectra 
in each globular cluster coadd is approximately equiva- 
lent to the number of spectra included in each ECHOS 
coadd. 

From the n MPMSTO spectra that remain after the 
application of our cuts in equatorial coordinate, radial 
velocity, and g — r color, we select with replacement 
m = 1,2,3, ...,ri spectra from the available data. We 
coadd the spectra and use the SSPP to derive [Fe/H] and 
[a/Fe] for that bootstrap coadded spectrum, and save 
the result. We repeat that process 100 times for each of 
m = 1, 2, 3, . . . , n to build up the distribution of SSPP 
[Fe/H] and [a/Fe] estimates for both globular clusters as 
a function of S/N. For M 13, the SSPP produces an es- 
timate of [Fe/H] = -1.7 ± 0.15 and [a/Fe] =0.3 ±0.15 
High-r esolu tion measurements of M 13 by iKraft et al.l 
(fT997l) and lCohen fc: Melended (l2005l) yielded [Fe/H] = 
-1.59 with [a/Fe] = 0.22 and [Fe/H] = -1.55 with 
[a/Fe] = 0.26 respectively. Likewise, for M 15 we find 
that [Fe/H] = -2.4 ± 0.2 and [a/F e] = 0.3 ± 0.2 High - 
reso lution measurements o f M 15 bv lSneden et al.l (|1997[ ) 
and iSneden eFall (|2000bD yielded [Fe/H] = -2.19 with 
[a/Fe] = 0.38 and [Fe/H] = -2.28 with [a/Fe] = 0.40 
respectively. We present the results for M 13 in Figure [3] 
and the results for M 15 in Figure U) This implies that at 
the resolution of the SEGUE spectra and for the range 
of Tcff' and \ogg included in the coadds, S/N is the dom- 
inant contribution to the MSE in estimating the average 
metallicity of a population using the ensemble of boot- 
strap coadds. At the same tim e, these results co mpare 
favor ably to tho se presented in lLee et al.l (|2008al lbl') and 
AUen de Prieto et al. (2008), as we measure the perfor- 
mance of the SSPP on a subset of the SEGUE data in a 
narrow range in g — r color. 

We summarize the precision and accuracy of our SSPP 
analysis of coadded spectra in Figure [5l The precision 
and accuracy of our SSPP analysis is a function of metal- 
licity; consequently, our estimates are most precise and 
accurate (~ 0.1 dex in [Fe/H] and [a/Fe]) for the most 
metal-rich populations, and least precise and accurate for 
the most metal-poor populations (~ 0.2 dex in [Fe/H] 
and [a/Fe]). Moreover, we find that at equivalent metal- 
licity, the MSE we compute for M 13 and M 15 based on 
coadd spectra created by coadding MPMSTO spectra in 
a finite range of effective temperature and surface grav- 
ity are in good agreement with the MSE computed in the 
ideal case of constant effective temperature and surface 
gravity for coadds based on noise-degraded MPMSTO 
spectra. As a result, we use the precision we derived for 
the SSPP metallicity analysis of the noise-added spectra 
to characterize the precision of our ECHOS metallicity 
measurements. For that reason, including both statisti- 
cal and systematic effects, our SSPP analysis of coadded 
SEGUE MPMSTO spectra produces estimates that are 
sufficiently precise and accurate to identify chemical dif- 
ferences between ECHOS and the kinematically smooth 
inner halo MPMSTO population on the order ~ 0.2 dex 
in both [Fe/H] and [a/Fe]. We use the expected MSE 
as a function of metallicity and S/N given in Figure [5] 
to determine whether an ECHOS is chemically distinct 
from the kinematically smooth halo MPMSTO popula- 
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tion along the same line of sight. 

3.5. The Metallicity of ECHOS 

For each of the three classes of ECHOS from lSOOl we 
consider for coaddition those spectra that correspond to 
stars within a radial velocity overdensity and ther efore 
consistent with ECHOS membership. As defined in lSOQl 
a radial velocity overdensity is a group of MPMSTO stars 
observed with radial velocities within a narrow range 
such that the group is extraordinarily unlikely to be ob- 
served if the underlying population were kinematically 
smooth. For the bin detections, we coadd the spectra of 
stars in the 20 km s~^ bin that contains the significant 
detection. For the peak detections, we coadd the spectra 
of stars within the measured width of the peak i n th e 
cumulative distribution function (the O statistic of lS09D . 
The width of the peak in 9 is the same as the velocity 
dispersion quoted for each ECHOS in Tables [2] and [S] 

As before, we select those stars within the 500 K wide 
bin in effective temperature that produces the highest 
S/N in the resultant coadded spectrum. We shift each 
candidate spectrum to = km s~^ and interpolate the 
spectrum and its inverse variance on to a common wave- 
length grid, and rescale each spectrum and its inverse 
variance to a common scale. From the n radial velocity 
selected, temperature controlled, radial velocity zeroed, 
interpolated, and scaled MPMSTO spectra, we sample n 
with replacement and coadd the spectra to create a single 
bootstrap realization of the resultant coadded spectrum 
distribution. For each ECHOS, we repeat the resampling 
process 100 times to create 100 bootstrap coadded spec- 
tra. We then run the SSPP on all of the bootstrap coad- 
ded spectra to determine the distribution of measured 
[Fe/H] and [a/Fe] for the ensemble of realizations. 

As a control, along the same line of sight as each 
ECHOS we perform the same steps on all the MPM- 
STO spectra that are associated with the kinematically 
smooth halo population. As a result, the stars in the 
control sample are in the same volume, were observed at 
the same time, and have a similar range of magnitude 
and S/N as the stars in the ECHOS. In that way, we 
can characterize the mean of [Fe/H] and [a/Fe] and esti- 
mate the error in our estimate of the mean for both the 
ECHOS and the kinematically smooth population along 
the same line of sight. With that information, any ob- 
served difference in composition is unlikely to result from 
systematic effects. Consequently, any observed chemical 
offset between an ECHOS and the kinematically smooth 
inner halo MPMSTO population along the same line of 
sight likely results from genuine chemical differences. 

We report the results of these calculations for all of 
the ECHOS discovered in S09 in Figure [6] through Fig- 
ure [11] and in tabular form in Table [T] through Table [31 
In Figure [6] through Figure ITl] the precision of our SSPP 
analysis is best characterized by the scatter apparent 
in results of the bootstrap coaddition process for both 
ECHOS and the kinematically smooth halo. If the two 
clouds of points do not overlap, than the two popula- 
tions are chemically distinct. We quantify this scatter 
and report the result in Table [T] through Table [Sj In all 
cases where the quoted metallicity of the smooth com- 
ponent is more iron-rich than typically associated with 
the smooth inner halo, the reason is because the ECHOS 
dominates the MPMSTO population along that line of 



sight (see Figures 2 through 11 of IS09I ). As a result, it is 
difficult to identify a large enough sample of stars in the 
kinematically smooth inner halo MPMSTO population 
for the equivalent analysis, so the quoted metallicity of 
the smooth halo can be significantly influenced by the 
metallicity of the ECHOS. 

4. DISCUSSION 
4.1. The Kinematic and Chemical Properties of ECHOS 

We searched for correlat ions between the dynamical 
properties of ECHOS from IS09I and the chemical prop- 
erties determined in this analysis. For the rest of this 
section, we examine the prop ertie s of the 21 Class II 
peak ECHOS from Table 3 of fSOOl. That sample of 21 
ECHOS is both the largest sample and the most repre- 
sentative sample of the inner halo ECHOS population. 
In general, we find that ECHOS that are iron-rich also 
have large radial velocity dispersions and are the most 
prominent ECHOS in that they have high number den- 
sities and are fractionally the largest contributors to the 
MPMSTO population along the line of sight where they 
were discovered; we plot these relations in Figure 1121 
Though large velocity dispersions are found only for the 
most metal-rich ECHOS, small velocity dispersions are 
found at all metallicities. The pro minence of metal- 
rich substructures was predicted by iFont et all (|2006D . 
as metal-rich substructures are preferentially produced 
by the most luminous progenitors. We find no signifi- 
cant correlation between [Fe/H] or [a/Fe] and distance. 

We initially identified ECHOS because they are kine- 
matically distinct from the kinematically smooth inner 
halo MPMSTO population. As we showed in Figure [6] 
through Figure m] nearly all ECHOS are also chemically 
distinct from the background smooth inner halo MPM- 
STO population along the same line of sight. We summa- 
rize this chemical distinctiveness in Figure 1131 As a pop- 
ulation, ECHOS are more iron-rich but less a-enhanced 
than the kinematically smooth background inner halo 
MPMSTO population. We showed in S09 that 10% of 
the inner halo (by volume) has 30% of its MPMSTO pop- 
ulation in ECHOS. Combined with the observation that 
ECHOS are metal-rich, these facts suggest that the most 
likely metallicity for an accreted star in the inner halo is 

[Fc/H] 1. At ECHOS [Fe/H] > -1, the apparent 

correlation between the iron metallicity of ECHOS and 
the iron metallicity of the smooth component of the halo 
results from the fact that iron-rich ECHOS are also the 
most prominent ECHOS (see Figure [T^ . That is, since 
the iron-rich ECHOS make up such a large fraction of the 
halo MPMSTO population along the line of sight where 
they were discovered, it is difficult to isolate a sample of 
MPMSTO stars in the kinematically smooth halo popu- 
lation. 

To asses the statistical significance of the observation 
that the population of ECHOS is more iron-rich but less 
a-enhanced than the kinematically smooth background 
inner halo MPMSTO population, we used a Monte Carlo 
simulation. Imagine that ECHOS and the kinematically 
smooth component of the halo really did have the same 
chemical composition. In that case, in Figure [T3I the de- 
parture of the points from the line y = x must result 
from imperfect observation, characterized by the error 
bars in the plot. Under that null hypothesis, we sam- 
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pie each ECHOS composition from a normal distribu- 
tion centered on the line y = x with standard deviation 
equal to the error in our measurement of the composition 
of the ECHOS. Likewise, we sample the composition of 
the smooth component along each line of sight from a 
normal distribution centered on the same point on the 
line y = X with standard deviation equal to the error in 
our measurement of the composition of the smooth com- 
ponent. We compute the signed, cumulative Euclidean 
distance of the entire population from the line y = x, 
and save the result. We repeat this process 10^ times. 
We find that under the null hypothesis, in no instance 
does the Monte Carlo simulation produce a cumulative 
distance of each point from the line y = x equal to the 
cumulative distance observed in the ECHOS [Fe/H] dis- 
tribution. Therefore the probability that the population 
of ECHOS has the same [Fe/H] as the smooth popula- 
tion along each line of sight in which we discovered an 
ECHOS is less than 10~^. Our [a/Fe] estimates are much 
less precise than our [Fe/H] estimates; nevertheless, the 
probability that the population of ECHOS has the same 
[a/Fe] as the smooth population along each line of sight 
in which we discovered an ECHOS is about 10^'^. The 
fact that ECHOS, as a population, are so chemically dis- 
tinct from the smooth background inner halo MPMSTO 
population along the same line of sight stron gly su pports 
the kinematic substructure identifications in 'SOS. 

The stars in ECHOS preferentially have apparent mag- 
nitudes that place them in the most dista nt half of our 
MPMSTO sample (jSchlaufman et al.ll2009l ). As a result, 
it is possible that the chemical distinctiveness of ECHOS 
is the result of a metallicity gradient in the inner halo. 
To test this hypothesis, we considered only those lines of 
sight that were not targeted at a known element of sub- 
structure and for which we had no significant ECHOS 
detection of any kind. In other words, those lines of 
sight which are dominated by the smooth component of 
the halo. For every MPMSTO star in that sample, we 
very roughly estimate the distance to each star accord- 
ing to the procedure described in Section 2 of S09. We 
then split the sample in half at about 14 kpc: the nearest 
50% of the MPMSTO population goes into the "close" 
subsample and the farthest 50% of the MPMSTO popu- 
lation goes into the "far" subsample. 

Then, for each subsample, we compute line of sight av- 
erage metallicities by the same spectral coaddition pro- 
cess described in Section 3. In that way, we end up with 
two estimates for the average chemical composition of 
MPMSTO stars in the smooth component along each 
line of sight where there is no significant substructure - 
an estimate for the "close" subsample and estimate for 
the "far" subsample. We find that the average metal- 
licity for the "close" subsample is [Fe/H] = —1.7 ±0.1 
while the average metallicity for the "far" subsample is 
[Fe/H] = —1.5 ±0.1. Meanwhile, the average metallicity 
of our ECHOS is [Fe/H] = -1.1 ± 0.1. For that rea- 
son, the apparent chemical distinctiveness of ECHOS is 
unlikely to be the result of a metallicity gradient in the 
inner halo. Indeed, though the halo is pot entially chem- 
ically inhomogeneous on large scales (e.g., iCarollo et al.l 
120071 [20Toh . our analysis is confined to a relatively small 
region in the inner halo. 

ECHOS are also chemically distinct from the other 



components of the Milky Way. The average level of star 
formation appears to have been mo re or less continu- 
ous in the thin disk for many Gyr ()Rocha-Pinto et all 
2000a|3), so thermonuclear SN have had plenty of 
time to contribute to the chemistry of the ISM. Con- 
sequently stars in the thin disk are typically close 
to solar in both [Fe/H] and [a/Fe]. The origin 
of thick disk stars is unclear, though they are uni- 
formly more g-enha nced than thin-disk s tars at con- 
stant [Fe/H] (e.g., iBensbv et al.l 120051 iReddv et all 
[200l IBensbv et ail 
our ECHOS in the 



20071 ). 



In Figure [14] we plot 
[a/Fe] plane along with 
[Fe/H] and [a/Fe ] estimates for individ ual st ar s from 
lEdvardsson et all (I1993D. [N isscn & Sc husteii (I1997D 



[Fe/H]- 



Hanson et al.' fl998'), 'Fulbright (2000), Prochask a et al 
(2000), Stcohcns&Bocsgaard (2002), Bcnsbv et al 
(2003), and Rcdd v et all (|200 3) and presented in 
Venn et al. (2004). We plot only those stars from 



Venn et aH (|2004[ ) that have a better than 50% assoca- 



tion probability with the thin disk, thick disk, or halo 
as indicated by their U, V, and W velocities. ECHOS 
are found in a region of the [Fe/H]- [a/Fe] plane sparsely 
populated-but not completely devoid of-^individual stars. 
In general, ECHOS are: (1) more iron-rich and less a- 
enhanced than halo stars, (2) more iron-poor than typ- 
ical thick-disk stars, and (3) more iron-poor and more 
a-enhanced than typical thin-disk stars. As a result, 
most stars associated with recent accretion events are 
more metal-rich than the average metallicity of the inner 
halo. In Figureim we plot the location of ECHOS in the 
[Fe/H]-[a/Fe] plane along with the location of individual 
giant stars in eight dSph galaxies reported in lKirbv et al.l 
poTol ). If ECHOS are the debris of a tidally disrupted 
dSph, Figure [15] indicates that the progenitor may have 
been comparable to the Sculptor or Leo I dS ph galaxies. 

Inte restingly, there may be a hint in the iVenn et al.l 
(|2004f ) compilation that those stars that are on retro- 
grade orbits are the halo stars with [a/Fe] most like the 
stars in the ECHOS (though the ECHOS have higher iron 
metallicity). Several studies using nearby halo star sam- 
ples have found correlations between increased scatter or 
peculiar elemental abundan ce patterns and e xtreme kine- 
matics (iCarney et al 1119971: iF ulbright 20 0^ llvans et al.l 
[200l IRoedered 12001 ICarollo et al.i ,20071. 12010^ . These 
studies find that metal-poor stars belonging to the dis- 
tant, outer halo originated in a more varied and/or inho- 
mogeneous environment, in support of the idea that the 
outer halo is assembled by more recent accretion of many 
low-mass systems. However, those studies are limited to 
[Fe/H] < — 1, and our results demonstrate that for stars 
currently in the inner halo region of the Galaxy, most 
stars accreted in the last 5 Gyr are not metal-poor, but 
instead have [Fe/H] ^ —1. Any comprehensive study 
of the halo accretion history must include stars over a 
broad range in metallicity. A direct comparison between 
our results and the metal-poor samples would require an 
analysis of selection biases in those samples that is be- 
yond the scope of this paper (for a discussion of the issues 
see iRoed erer 2009), but we note that the sample of halo 
stars in which the ECHOS were identified was selected 
using ultraviolet excess, which wo uld te nd to bias our 
sample against metal-rich stars (see IS09I for details) . 

In Figure I16| we compare the velocity dispersions of 
ECHOS from lS09l with the median radial velocity pre- 
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cision of the MPMSTO stars in each ECHOS. In this 
case, th e rad ial velocity precision as a function of S/N we 
used in IS09I assumed a population metallicity of [Fe/H] 
~ —1.6. SEGUE radial velocity precision becomes bet- 
ter at higher metallicities, so the fact that ECHOS have 
[Fe/H] ^ —1.0 may indicate the actual radial velocity 
precision is somewhat better. In any case, we find many 
of our ECHOS have velocity dispersions that are substan- 
tially larger than the expected radial velocity precision. 

4.2. The Origin of ECHOS 

The stars in ECHOS are more iron-rich and less a- 
enhanced then the typical inner halo MPMSTO star pop- 
ulation. They are also more iron-poor than typical thick- 
disk stars but more a-enhanced than typical thin-disk 
stars. The high [Fe/H] metaUicity of ECHOS almost 
certainly rules out ultrafaint dSph galaxies as ECHOS 
progenitors, as ult rafaint dSph galax ies have average 
[Fe/H] < -2 fe.g.. IKirbv eTalll2008bD . If ECHOS are 
the tidal debris of one or more dSph galaxies, the dSph 
luminosity-metallicity relation (e.g.. IKirbv et al.]l2008b( ) 
implies a progenitor luminosity of L ~ 10^ Lq to produce 
a characteristic iron metallicity [Fe/H] ^ —1.0. That lu- 
minosity combined with a reasonable dSph stellar mass 
to light ratio Mtf/L y ^ 3 and t otal mass to light ratio 
Mtot/Lv ~ 10 (e.g.. lMateolll998| ) implies the accretion of 
a progenitor with stellar mass M,, ~ 3 x 10^ M© and total 
mass Mtot ^ 10^ Af©, comparable to Local Group mem- 
bers NGC 147 or NGC 185. Radial velocity substruc- 
tures persist for many crossing times, and because the 
typical crossing time in the inner halo is 50 Myr ECHOS 
shou ld be observable for up to 5 Gyr (see Section 5.1 of 
IS09I ). As a result, if the debris of an ECHOS progeni- 
tor is still visible as cold radial velocity substructure, the 
accretion event must have occurred in the last ~ 5 Gyr, 
or equivalently, since z = 0.5. Such accretion events are 
common in the accretion histories of Milky Way analog 
halos, and a typical halo might have experienced 10 such 
accretion events si nce z = 0.5 (e.g., Stewart et al. 2008; 
iFakhouri e tliL 20101. Again, the distribution of ECHOS 
in the [Fe/H] [a/Fc] plane is similar to the distribution 
in [Fe/H]-[a/Fe] plane of individ ual giant stars in the 
Sculptor dSph (iKirbv et al.llMol) . 

The radial velocity dispersion distribution of ECHOS 
may also be a clue to their origin. Collisionless dynam- 
ics implies that as the position space component of the 
phase-space distribution of a tidally disrupted stellar sys- 
tem becomes hotter, the velocity component of its phase- 
space distribution must n ecessarily become colder (e.g., 
iBinnev fc Tremaind 119871 ). In other words, as a tidally 
disrupted stellar system spreads across the sky its veloc- 
ity dispersion measured over a small patch of sky must 
decrease with time. For that reason, the radial veloc- 
ity dispersion of an element of cold halo substructure 
is a lower limit on the radial velocity dispersion of its 
progenitor. At the same time, multiple wraps of stream- 
like substructure might give the appearance of a single 
substructure with large velocity dispers ion, though there 
is both observational e vidence (e.e.. iBell et al.l I2008D 
and theoretical modeling ([Johnston et al.ll2008D that sug- 
gests stream-like substructures are rare in inner halo. If 
ECHOS really are shell-like substructures on radial orbits 
as opposed to stream-like substructures on tangential or- 



bits as proposed bv lS09l . measurement over a small patch 
of sky might still intersect orbits in a range of orbital 
phase and therefore produce large velocity dispersions 
that are not representative of the velocity dispersion of 
the progenitor. We argue that this is unlikely, however, 
as that would imply that the radial velocity dispersion 
within a single ECHOS is a function of apparent mag- 
nitude, a tre nd no t observed in ECHOS (see Figures 2 
through 11 in lS09D . With those caveats in mind, typical 
classical dwarf spheroidal (dSph) galaxies have radial ve- 
locity dispersions ct '--^ 10 km s~^ (e.g.. iMat co 1998). In- 
deed, three of our ECHOS even have velocity dispersions 
comparable to that observed in the Small Mag ellanic 
Cloud (SMC) in which IHarris fc Zar"itsWl (l200fif l mea- 
sured a « 27.5 ± 0.5 km s~^. The median radial velocity 
dispersion of ECHOS ct ~ 20 km s~^ is also characteristic 
of a dSph like NGC 147 or NGC 185. This observation 
implies an even more massive progenitor though, as the 
fact that ECHOS are not associated with surface bright- 
ness substructure suggests that the velocity dispersion 
has had time to cool measurably. 

The observation that ECHOS are less a-enhanced than 
typical inner halo stars suggests that the star formation 
timescale in the ECHOS progenitor was long relative to 
the star formation timescale in t he massive progen itor of 
the bulk of the inner halo (e.g.. iMcWilliam et al.| [l995). 
so long as the stellar initial mass function is not too dif- 
ferent in the two environments. This relatively extended 
star formation timescale might disfavor globular clusters 
as ECHOS progenitors, though some globular clusters 
have [g/Fe] ^ 0.2 (e.g.. lCarnevl ll99l iBrodie fc Stradeii 
l2006tlCarretta et al.|[2C)09D . On the other hand, the char- 
acteristic ECHOS [Fe/H] 1.0 falls between the two 

peaks in the observed bimoda l Milky Way globular clus- 
ter [Fe/H] distribution (e.g., lArmandroff k ZinnI 119881 : 
IBrodie fc Straded 120061 ). Additionally, globular clusters 
that metal-rich ar e very rare in the inner halo (e.g., 
IGeisler et al.]l2007l) . In any case, the large radial velocity 
dispersion of ECHOS indicates that ECHOS are unlikely 
to be the debris of tidally disrupted globular clusters, 
as globul ar clusters have velocity dispersions cr ~ 5 km 
s~^ (e.g.. iMandushev et al.l[T99l[ ). Core-collapse globu- 
lar clusters have higher velocity dispersions, but are also 
resistant to disruption. 

One last possibility is that the stars in ECHOS formed 
in the nascent disk of the Milky Way and were scat- 
tered into the inner halo during an accretion event. Re- 
cent cosmological m odels of Milky Way formation by 
iZolotov et al.l ()2009( ) suggest that scattering into the in- 
ner halo does contribute to the stellar population in the 
inner halo. Moreover, the accretion of dSph galaxies 
more massive then those suggested by the characteris- 
tic ECHOS [Fe/H] ~ —1.0 can scatter sub stantial num- 
bers of disk stars into the inner halo (e.g., iPurcell et al.l 
l2010f ). However, models of the chemical evolution of the 
disk of the Milky Way suggest that the typical chemical 
composition of ECH OS is difficu lt to explain with a disk 
population (Schonr ich fc Binnev 2009a, b). In addition, 
IZolotov et~al.i (,20lO| ) recently suggested that disk stars 
scattered into the halo should be more a-enhanced than 
inner halo stars at constant [Fe/H]. Improved resolution 
and better treatments of star formation and feedback in 
theoretical models may yet determine whether scattering 
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of disk stars into the inner halo by low-mass accretion 
events is a plausible origin of ECHOS. 



5. CONCLUSION 

The metal-poor main sequence turnoff stars in the 
elements of cold halo substructure in the inner halo 
of the Milky Way found to be kinematically distinct 
fro m the kinematica l ly sm ooth inner halo population 
by iSchlaufman et al.l (|2009l ) are also chemically distinct 
from the smooth inner halo. As a population, ECHOS 
are more iron-rich and less a-enhanced than the MPM- 
STO stars in the kinematically smooth halo population 
along the same line of sight. ECHOS are chemically dis- 
tinct from both the thin and thick disk as well, as they are 
more iron-poor than average thick-disk stars and both 
more iron-poor and more a-enhanced than average thin- 
disk stars. Kinematically, the typical velocity dispersion 
of the ECHOS population is cr ~ 20 km s~^, though 
ECHOS have a range of radial velocity dispersions ex- 
tending from the floor of cr ~ 10 km s""'^ set by the pre- 
cision of SEGUE radial velocities at the apparent mag- 
nitude of the stars in ECHOS to the maximum observed 
value of (T ^ 30 km s^^. If ECHOS are the result of the 
tidal disruption of an accreted dwarf spheroidal galaxy or 
galaxies, the high iron metallicity [Fe/H] ~ —1.0, large 
velocity dispersion cr 20 km s~^, and the lack of cor- 
responding surface brightness substructure imply the ac- 
cretion of a Aftot ^ 10^ Mq halo sometime in the last 5 
Gyr since z — 0.5. In addition, the high iron metallicities, 
low a-enhancements, and large velocity dispersions are 
difficult to reconcile with globular clusters as the progen- 
itors of ECHOS. One final possible explanation for these 
observations is that the stars in ECHOS formed in the 
nascent disk of the Milky Way and were subsequently 
scattered into the inner halo by cosmologically common 
low-mass accretion events, though models of the chem- 
ical evolution of the Milky Way's disk cast some doubt 
on this scenario as well. 
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APPENDIX 

A. DETAILED DESCRIPTION OF COADDITION ALGORITHM 

The selection of the input sample of spectra to the coaddition process is different in each of the three contexts 
described in Sections 3.3, 3.4, and 3.5. In all cases, we select a bootstrap sample of spectra, coadd the spectra in the 
bootstrap sample into a single coadded spectrum, and then use the SSPP to analyze the coadded spectrum. That 
analysis produces a single estimate of the mean metallicity of the population of stars from which the bootstrap sample 
was drawn. We draw many bootstrap samples, derive many estimates of the mean metallicity of the population, and 
use that distribution to quantify both the mean metallicity of the population and our uncertainty in that estimate. 

In Section 3.3, we estimate the accuracy and precision of the SSPP on coadded spectra in a range of metallicity 
(-2.41 < [Fe/H] < -0.23) and S/N (10 < S/N < 100). For each of eight MPMSTO stars with metallicity known from 
high resolution spectroscopy, the input sample is a collection of 640 spectra produced by adding noise to the original 
high S/N SEGUE spectrum of that star. There are no radial velocity or effective temperature cuts applied in this 
case. 

In Section 3.4, we estimate the accuracy and precision of the SSPP on coadded spectra produced by bootstrap 
sampling from the MPMSTO populations of the globular clusters M 13 and M 15. In this case, though the individual 
spectra belong to stars with very similar metallicities, the spectra will correspond to stars with different radial velocities 
and effective temperature. As a result, we can use the spectra associated with M 13 and M 15 MPMSTO stars to 
quantify the effect of coadding stars in a finite range of effective temperature, surface gravity, and radial velocity. 
Specifically, we include in the sample of input spectra only those spectra corresponding to MPMSTO stars with 
equatorial coordinates within the projected tidal radius of each cluster, with radial velocities close to the systematic 
radial velocity of the cluster, and within the 500 K range of effective temperature that produces the highest S/N in 
the final coadded spectrum. 



10 



Schlaufman et al. 



In Section 3.5, we use our algorithm to estimate both the average metallicity of MPMSTO stars in ECHOS and 
in the kinematically smooth component of the halo along the same line of sight. We include in each ECHOS sample 
of input spectra only those spectra corresponding to MP MSTO stars with radial velocities within the radial velocity 
overdensity that defines the ECHOS (as defined in IS09I and reported in Tables [1] [21 and|3]). We include all other 
MPMSTO spectra obtained along the line of sight in the kinematically smooth halo input sample. In both cases, we 
use the same effective temperature criteria used in the globular cluster case. 

The following eight steps describe how a single bootstrap coadded spectrum can be produced from the input samples 
described above. 

1. Use resampling with replacement to select N spectra from the N available spectra to create a bootstrap sample 
of spectra to synthesize into a single coadd spectrum. 

2. Let fi,j{X'j) and gi.j(A^) denote the fiux and inverse variance from the SDSS spectroscopic reduction pipeline. 
The subscript i refers to the ith star in the coadd and the subscript j refers to the jth wavelength bin. 



3. Shift each spectrum to rcdshift z = such that 



A' 



A, = ^^ (Al) 

^ 1+z ^ ' 

and denote the radial- velocity zeroed spectrum and inverse variance fi.j{)^j) and 9i,j{^j) respectively. 

4. Use natural cubic spline interpolation to interpolate fi,j{Xj) — /i(A) and gi,ji^j) 5i(A). 

5. Sample the continuous functions fi{X) and 5i(A) on to a common grid in Afc for every star from Ai = 3850 A to 
A2 = 9000 A in n 0.5 A increments. Set any fi^k{\k) < to fi^ki^k) = and any gi^ki^k) < to gi^ki^k) — 10^^. 

6. To this point, each star has a different total flux in its spectrum, so rescale fi and gi to the same total fiux while 
conserving the S/N ratio in each bin. In other words, 

A,= M\)d\ (A2) 

/»,fc(Afc) = ^Afe(Afe) (A3) 

g^^k{\k)^Alg^,k{\k)■ (A4) 

7. Loop over all n wavelength bins to create the coadd spectrum and its inverse variance using the weighted-mean 
such that 

Wi = {5ii,fc(A/c),52,fc(Afe),. . . ,5w,/c(Afe)} (A5) 

-f (\ \ UJjfi^ki^k) 

Jk[^k)^ — jv (Ab) 

N 

9ki^k)=J2'^^- (A^) 

i=l 

8. The coadded spectrum /^(Afc) and its inverse variance gj^{Xk) are then processed by the SSPP to produce a 
single estimate of the mean metallicity of the input sample. 

A single iteration of these eight steps produces a single estimate of the mean metallicity of the input sample, 
corresponding to a single point in Figures [T] through 2] and Figures [B] through 1111 We repeat these eight steps many 
times to produce many estimates of the mean metallicity of the input sample. We use that distribution to estimate 
both the true mean metallicity of the input sample and our uncertainty in that estimate. 
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Figure 1. Ensemble of bootstrap coadd realizations of single MPMSTO stars giving SSPP [Fe/H] precision and accuracy as a function of 
S/N for that s tar. The [Fe/H] and T^g estimates given for each star were derived from Hobby-Eberly Telescope High Resolution Spectrograph 
observations lAllende Prieto et al.,,.200 8'1: we also include the SSPP [Fe/H] estimate for the original, single high S/N spectrum. Each point 
in the plot corresponds to the result of a SSPP analysis of a single bootstrap coadded spectrum made up of many model low S/N MPMSTO 
star spectra created by degrading a high S/N SEGUE MPMSTO star spectrum with a detailed noise model. The mean square error (MSE 
= bias^ + variance) of our SSPP analysis of codded spectra ranges from 0.05 dex in [Fe/H] for the most iron-rich stars to 0.20 dex in 
[Fe/H] for the most iron-poor stars. 
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Figure 2. Ensemble of bootstrap coadd realizations of single MPMSTO stars giving SSPP [a/Fe] precision and accuracy as a function of 
S/N for that star. The panels are in the same order as in Figure [T] The [a/Fe] and T^g estimates given for each star were derived from 
Hobby-Eberly Telescope High Resolution Spectrograph observations IIAUende Prieto et al.ir2008 ^: we also include the SSPP [o/Fe] estimate 
for the original, single high S/N spectrum. Each point in the plot corresponds to the result of a SSPP analysis of a single bootstrap coadded 
spectrum made up of many model low S/N MPMSTO star spectra created by degrading a high S/N SEGUE MPMSTO star spectrum 
with a detailed noise model. The MSE of our SSPP analysis of codded spectra ranges from 0.05 dex in [a/Fe] for the most iron-rich stars 
to 0.20 dex in [a/Fe] for the most iron-poor stars. 
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Figure 3. Ensemble of bootstrap coadd realizations of single MPMSTO stars giving SSPP [Fe/H] and [a/Fe] precision and accuracy as a 
function of S/N for M 13. Each point in the plot corresponds to a single bootstrap coaddition of SEGUE MPMSTO spectra corresponding 
to stars with equatorial coordinates within the projected tidal radius of M 13, with radial velocities within 15 km s~^ of the systematic 
radial velocity of the cluster, and with g — r colors in an interval that corresponds to a range in effective temperature as noted in each panel. 
Note that the precision and accuracy of the SSPP analysis of the coadded spectra is not a strong function of the g — r color range (and 
therefore T^g range) of the MPMSTO stars included in the coadd. Combined with the fact that coadding in 500 K bins always leads to a 
higher S/N coadded spectrum, we coadd spectra corresponding to MPMSTO stars in the bin of 3 — r color corresponding to the T^g bin of 
width 500 K centered on the T^g value that gives the highest S/N in the resultant coadded spectrum. In this way, the SSPP produces an 
estimate for the metallicity of M 13 of [Fe/H] = — 1 .7±0.15 and [a/Fe] = 0.3±0.15 in rough agreement with previous analyses iKraft et al.l 
\T99% ICohen fc Melende2|[2005l: IKirbv et al]|2008al') . 
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Figure 4. Ensemble of bootstrap coadd realizations of single MPMSTO stars giving SSPP [Fe/H] and [a/Fe] precision and accuracy as a 
function of S/N for M 15. Each point in the plot corresponds to a single bootstrap coaddition of SEGUE MPMSTO spectra corresponding 
to stars with equatorial coordinates within the projected tidal radius of M 15, with radial velocities within 25 km s~^ of the systematic 
radial velocity of the cluster, and with g — r colors in an interval that corresponds to a range in effective temperature as noted in each 
panel. In this way, the SSPP pro duces an estimate for the m etallicity of M 15 o f [Fe/H] = —2.4 ± 0.2 and [a/Fe] = 0.3 it 0.2 in rough 
agreement with previous analyses HSneden et al.|[T997l . [2000al lb HKirbv et al.|[2008al 'l . 
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Figure 5. The mean square error (MSE = bias^ + variance) of our SSPP analysis of coadded MPMSTO spectra as a function of S/N 
and metallicity. The blue curves correspond to our analysis of the noise-added spectra, while the green curves correspond to our analysis 
for M 13 and M 15. The widths of the curves are common between both panels in that the curve for a single star is plotted with the 
same width in both panels. In both panels, the most solar-like compositions are always the lightest shades of blue. Note that the SSPP 
produces more precise [Fe/H] and [o/Fe] estimates at common S/N for the most metal-rich MPMSTO stars. In that case, the abundance 
of metal lines in the moderate resolution SEGUE spectra of metal-rich stars permits precise [Fe/H] estimates. Also note that the MSE of 
our analysis is more strongly affected by population metallicity than it is affected by the coaddition of MPMSTO star spectra with a small 
but finite range of T^f{ and log g. For that reason, the expected MSE of our analysis is best determined by our analysis of the noise-added 
spectra. Ultimately, including both statistical and systematic effects, our SSPP analysis of coadded SEGUE MPMSTO spectra produces 
estimates that are precise and accurate enough to identify differences on the order of 0.2 dex in both [Fe/H] and [o/Fe]. Left: SSPP [Fe/H] 
MSE as a function of S/N and metallicity. Right : SSPP [a/Fe] MSE as a function of S/N and metallicity. 



18 



Schlaufman et al. 



n — ' — ' — ' — ' — r 


— 1 — 1 — 1 — r 
























ECHOa <|FbIH1» = -1.de _ 



Figure 6. Ensemble of bootstrap coadd realizations giving SSPP [Fe/H] estimates for ECHOS from 'S09I using bootstrap coaddition of 
individual SEGUE MPMSTO star spectra. Each blue point corresponds to a SSPP [Fe/H] estimate for a single ECHOS bootstrap coadd, 
while the dark gray points are SSPP [Fe/H] estimates from an equivalent bootstrap coaddition of MPMSTO star spectra with radial 
velocities indicating membership in the kinematically smooth inner halo MPMSTO population along the same SEGUE line of sight that 
hosts the ECHOS. MPMSTO stars correspond to an ECHOS if they have a radial velocity consistent with radial velocity overdensities as 
determined by S09; we identify all other stars along the same line of sight as the kinematically smooth inner halo MPMSTO population. 
Of those MPMSTO stars that are radial velocity members of an ECHOS, we include in the coadd only those spectra that correspond to 
stars that have g — r colors corresponding to tlie Tj.ff bin of width 500 K centered on the T^ff value that gives the highest S/N in the 
resultant coadded spectrum. In the same way, we choose an optimal T^ff bin of width 500 K from which we select spectra for coaddition 
from the kinematically smooth inner halo MPMSTO population. The light gray polygon is centered on the mean metallicity of the smooth 
component of the inner halo along each SEGUE line of sight and indicates the expected MSE of the SSPP at the metallicity of the smooth 
component as a function of S/N from our bootstrap coaddition of high S/N spectra degraded with the noise model (see Figure^. Therefore, 
if the cloud of blue points does not coincide with the gray polygon, then the ECHOS is chemically distinct from the kinematically smooth 
population along that line of sight. In all cases where the quoted metallicity of the smooth component is more iron-rich than typically 
associated with the smooth inner halo, the reason is because the ECHOS dominates the MPMSTO population along that line of sight (see 
Figures 2 through 11 of lSOg) . Left: SSPP [Fe/H] analysis for ECHOS B-1, B-2, B-3 (Monoceros), B-4 tCrillmair fc Dionatosii2006 ). and 
B-5 from Tablejl] Right: SSPP [Fe/H] analysis for ECHOS B-6, B-7, and B-8 from Table[T]as well as PCI-1 and PCI-2 from TableEl 
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Figure 7. Left: SSPP [Fc/H] analysis for ECHOS PCI-3, PCI-4 (Monoceros), PCI-5 HGrillmair &: Dionato£|[200g ) . PCI-6, and PCI-7 from 
TabloE] Right: SSPP [Fo/H] analysis for ECHOS PCI-8 and PCI-9 (Monoceros) from Table|2]as well as PCII-2, PCII-4, and PCII-5 from 
Table [S] See the caption to Figure |6] for a detailed description of this type of figure. 
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Figure 8. Left: SSPP [Fe/H] analysis for ECHOS PCII-6, PCII-7, PCII-9, PCII-13, and PCII-14 from Table [3] Right: SSPP [Fe/H] 
analysis for ECHOS PCII-15, PCII-16, PCII-17, PCII-18, PCII-19 from Tabled See the caption to Figure E] for a detailed description of 
this type of figure. 
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Figure 9. SSPP [a/Fc] estimates for ECHOS from'SO^ using bootstrap eoaddition of SEGUE MPMSTO star spectra. Left: SSPP [a/Fc] 
analysis for ECHOS B-1, B-2, B-3 (Monoceros), B-4 (Grillmair & Dionatos 2006), and B-5 from Table[T] Right: SSPP [a/Fe] analysis for 
ECHOS B-6, B-7, and B-8 from Table[T]as well as PCI-1 and PCI-2 from Tablc[2| See the caption to Figure[6]for a detailed description of 
this type of figure. 
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Figure 10. Left: SSPP [a/Fe] analysis for ECHOS PCI-3, PCI-4 (Monoceros), PCI-5 HGrillmair &: DionatosI [2009 ) . PCI-6, and PCI-7 
from TableE] Right: SSPP [a/Fe] analysis for ECHOS PCI-8 and PCI-9 (Monoceros) from Table[2]as well as PCH-2, PCn-4, and PCII-5 
from Table |3] See the caption to Figure [6] for a detailed description of this type of figure. 
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Figure 11. Left: SSPP [a/Fe] analysis for ECHOS PCII-6, PCII-7, PCII-9, PCII-13, and PCII-14 from Table E] Right: SSPP [a/I 
analysis for ECHOS PCII-15, PCII-16, PCII-17, PCII-18, PCII-19 from Table El See the caption to Figure [6] for a detailed description 
this type of figure. 



24 



Schlaufman et al. 



Ui 

E 

ID 



I — • — I 
I * » I 
.1 • 



I — • — I 



,LmJ I 



-2.0 -1.5 

ECHOS [Fe/H] 



n — — ' — ' — ' — — 

~# ECHOS 


T 1 1 1 \ \ 1 1 1 \ \ \ p 


- 


1 ■ 1 






I-«H 


' — • — \ lafeir^ 











-2.0 -1.5 

ECHOS [Fe/H] 



Figure 12. Kinematic properties of ECHOS as a function of metallicity. Left: Iron-rich ECHOS can have large v elocity disp ersions. 
Right: Iron-rich ECHOS have the highest space densities; this prominence of metal-rich substructures was predicted in lFont et all {2006). 
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Figure 13. Comparison of tlie metallicity of ECHOS to tiie metallicity of ttie MPMSTO population in the kinematically smooth component 
of the inner halo along the same line of sight where the ECHOS was discovered. Left: ECHOS are more iron-rich than the MPMSTO 
stars associated with the kinematically smooth inner halo. Moreover, ECHOS are so iron-rich that an origin within one or more now 
tidally disrupted ultrafaint dwarf galaxies is unlikely. That is, the mass— luminosity relation of Milky Way dwarf spheroidal galaxies (e.g., 
IKirbv et al.l l2008H) requires a progenitor luminosity L ~ 10® Lq to reach a mean iron metallicity [Fe/H] ~ —1.0. Lines of sight hosting an 
ECHOS for which we infer a smooth component metallicity [Fe/H] > — 1.6 are those lines of sight along which the ratio of the number of 
MPMSTO stars with radial velocities consistent with ECHOS membership Us to the total number of MPMSTO stars with radial velocity 
measurements Ns is non-negligible, or Us/Ns > 0. In other words, lines of sight where ECHOS are a significant contributor to the MPMSTO 
population tend to have inferred smooth component metallicities [Fe/H] > —1.6. As a result, those apparently enhanced iron metallicities 
in the smooth component are likely due to the presence of stars that are outside the radial velocity overdensity that defines the ECHOS, 
but that are still associated with ECHOS. Right: ECHOS are less «-enhanced than the MPMSTO stars associated with the kinematically 
smooth inner halo, so (modulo any changes in the initial mass function) the MPMSTO stars in ECHOS formed in environments in which 
the star formation tim escale was long relative to the star formation timescale in the massive progenitor of the bulk of the inner halo (e.g., 
[Robertson et al.|[2005tl . 
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Figure 14 . Distribution of ECHOS in the [Fe/H]-[a/Fe] plane along with individual stars collected from the literature as presented in 
IVenn et al.l lj2004). The ECHOS mostly fall in a region of that plane that is relatively sparsely occupied by-but not completely barren 
of-individual stars. In general, the ECHOS are more iron-poor than thick-disk stars and more iron-poor and o-enhanced than typical 
thin-disk stars. 
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Figure 15. Distribution of ECHOS in the [Fe/H]-[o/Fe] plane along with individual stars in dSph galaxies from IKirbv et aP l|20igV In 
each panel, the gray points in the background represent [Fe/H] and [a/Fe] measurements for individual dSph stars. The blue points in 
the foreground are our [Fe/H] and [a/Fe] measurements for individual ECHOS. Left to right and top to bottom, the dSph galaxies are: 
Sculptor, Fornax, Leo I, Sextans, Leo II, Canes Venatici I, Ursa Minor, and Draco. ECHOS are plausibly associated with a progenitor 
comparable to Sculptor or Leo I. 
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Figure 16. Velocity dispersion of ECHOS as a function of the median radial velocity precision of the stars in the ECHOS. The fact 
that the points are above the line y = x suggests that the observed velocity dispersion is not only a function of the limited radial velocity 
precision for faint MPMSTO stars. The velocity dispersion observed in an ECHOS is a lower limit on the velocity dispersion of any 
gravitationally bound progenitor. This results from the fact that the phase-space distribution of a disrupted stellar system becomes colder 
with time, combined with the fact that each SEGUE line of sight would only intersect a small fraction of an or bit of an elem ent of halo 
substructure. Typical classical dwarf spheroidal galaxies have radial velocity dispersions ~ 10 km (e.g.. IMateol ll998). while the 
SMC has a velocity dispersion of cr Ri 27.5 ± 0.5 km s"'^ fe.e.. [Harris fc Za.ritskvij2006i) . The origin of ECHOS through the gravitational 
disruption of a progenitor with a velocity dispersion cr > 20 km s"'^ implies a progenitor mass Aftot ^ 10^ Mq. 
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Table 1 

Properties of ECHOS 



ID 


RA 


Dec 


1 


b 


bplatc 


fplate 






Vol 


d 


err 






[Fc/H]b 


err 


[Fe/H]s 


err [o 




(dcg) 


(dcg) 


(dcg) 


(dcg) 










(kpc3) 


(kpc) 


(kpc) 


(km s-i) 




(dex) 


(dcx) 


(dex) 


(dex) 1 


B-1 


21.3 


39.6 


130.0 


-22.8 


2043 


2063 


34 


228 


2.20 


18.4 


3.2 


-130 


12 


-0.70 


0.12 


-0.99 


0.18 


B-2 


17.0 


0.0 


132.0 


-62.6 


2313 


2328 


109 


561 


3.48 


10.8 


6.6 


-170 


20 


-1.33 


0.14 


-1.71 


0.25 


B-3'' 


39.7 


28.2 


150.0 


-29.0 


2442 


2444 


59 


265 


2.30 


11.2 


2.0 


-50 


17 


-0.87 


0.09 


-0.90 


0.15 


B-4b 


163.8 


48.0 


162.4 


59.2 


2390 


2410 


150 


672 


3.86 


6.6 


2.8 


-130 


25 


-2.12 


0.11 


-1.82 


0.13 


B-5 


129.6 


53.9 


164.3 


37.2 


2316 


2331 


93 


425 


3.33 


9.6 


2.5 


-10 


20 


-1.02 


0.09 


-1.51 


0.17 


B-6 


124.5 


38.0 


183.4 


32.6 


2670 


2674 


83 


514 


3.46 


9.9 


1.1 


30 


17 


-0.97 


0.06 


-1.10 


0.13 


B-7 


132.6 


6.1 


221.5 


29.2 


2317 


2332 


69 


470 


3.04 


11.0 


1.3 


70 


17 


-0.93 


0.10 


-1.50 


0.18 


B-S'' 


134.0 


3.2 


225.2 


29.0 


2888 


2913 


74 


514 


3.50 


11.0 


1.7 


90 


19 


-1.08 


0.05 


-1.33 


0.14 



Note. — Kinematic and metallicity data for all bin detection ECHOS from S09|. The columns are; right ascension, declination, galactic longitude, 
plate number, SEGUE faint plate number, number of SEGUE MPMSTO spectra along the indicated line of sight, number of photometric MPMSTC 
line of sight, volume scanned along the indicated line of sight, median distance of the stars in the ECHOS, the error in that estimate, the radial velocil 
MPMSTO stars kinematically associated with the ECHOS, the iron metallicity of the ECHOS, the error in that estimate, the iron metallicity of the 1 
of the halo along the indicated line of sight, the error in that estimate, the a-enhancement of the ECHOS, the error in that estimate, the «-enhancem 
component of the halo along the indicated line of sight, and the error in that estimate. 

^ Monoceros Stream 

[Grillmair Sz DionatosI | |2003) Stream 
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Table 2 

Properties of ECHOS 



ID 


RA 


Dec 


1 


b 


bplate 


fplate 


Ns 


Np 


Vol 


d 


err 




cr 


err 


ns 


[Fc/U]e 


e: 




(deg) 


(deg) 


(deg) 


(deg) 










{kpc3) 


(kpc) 


(kpc) 


(km s-i) 


(km s^^) 


(km s-i) 




(dex) 


(d, 


PCI-1 


214.8 


56.4 


100.7 


56.8 


2447 


2462 


122 


673 


2.79 


8.7 


2.3 


-328 


15.1 


11.5 


8 


-1.73 





PCI-2 


20.0 


31.7 


130.0 


-30.8 


2041 


2061 


93 


349 


2.60 


13.4 


3.3 


-125 


22.0 


17.8 


20 


-0.95 





PCI-3 


21.3 


39.6 


130.0 


-22.8 


2043 


2063 


34 


228 


2.20 


18.7 


3.5 


-121 


10.5 


19.4 


13 


-0.65 





PCI-4'' 


39.7 


28.2 


150.0 


-29.0 


2442 


2444 


59 


265 


2.30 


11.4 


2.6 


-57 


10.2 


12.8 


22 


-0.87 





PCI-5'' 


163.8 


48.0 


162.4 


59.2 


2930 


2410 


150 


672 


3.86 


6.6 


2.5 


-132 


11.7 


4.6 


22 


-2.06 





PCI-6 


129.6 


53.9 


164.3 


37.2 


2316 


2331 


93 


425 


3.33 


10.0 


2.4 


-13 


10.9 


9.8 


20 


-1.03 





PCI-7 


124.5 


38.0 


183.4 


32.6 


2670 


2674 


83 


514 


3.46 


10.1 


1.5 


29 


19.0 


10.6 


16 


-0.91 





PCI-8 


132.6 


6.1 


221.5 


29.2 


2317 


2332 


69 


470 


3.04 


10.9 


1.6 


71 


13.2 


11.4 


15 


-0.95 





PCI-9'' 


134.0 


3.2 


225.2 


29.0 


2888 


2913 


74 


514 


3.50 


10.3 


1.4 


85 


14.9 


4.0 


17 


-1.01 





Note. 


— Kinematic and metallieity data for all elass I peak detection ECHOS from [50S|. The columns 


are: right 


ascension, declination, galactic 



SEGUE faint plate number, number of SEGUE MPMSTO spectra along the indicated line of sight, number of photometric MPMSTO candidates aL 
indicated line of sight, median distance of the stars in the ECHOS, the error in that estimate, the radial velocity of the ECHOS, the velocity dispersi 
stars in the ECHOS, the number of MPMSTO stars kinematically associated with the ECHOS, the iron metallieity of the ECHOS, the error in that 
component of the halo along the indicated line of sight, the error in that estimate, the o-enhancement of the ECHOS, the error in that estimate, the t 
the halo along the indicated line of sight, and the error in that estimate. 

Monoceros Stream 

^ IGrillmair fc DionatosI 1 12003) Stream 
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Table 3 

Properties of ECHOS 



ID 


RA 


Dec 


1 


b 


bplate 


fplate 






Vol 


d 


err 






err 


ns 


[Fe/H]B 




(dcg) 


(deg) 


(deg) 


(deg) 










(kpc3) 


(kpc) 


(kpc) 


(km s~^) (km 


s-i) 


(km s^-*-) 




(dcx) 


PCII-1 


20.0 


31.7 


130.0 


-30.8 


2041 


2061 


93 


349 


2.60 


13.4 


3.3 


-125 


22.0 


17.8 


20 


-0.92 


PCII-2 


20.0 


31.7 


130.0 


-30.8 


2041 


2061 


93 


349 


2.60 


12.7 


3.4 


-98 


24.8 


16.2 


18 


-0.99 


PCII-3 


21.3 


39.6 


130.0 


-22.8 


2043 


2063 


34 


228 


2.20 


18.7 


3.5 


-121 


10.5 


19.4 


13 


-0.67 


PCII-4 


91.8 


83.5 


130.0 


25.7 


2540 


2548 


47 


223 


2.45 


12.7 


1.9 


-95 


19.9 


12.0 


13 


-0.81 


PCII-5 


17.0 


0.0 


132.0 


-62.6 


2313 


2328 


109 


561 


3.48 


12.0 


6.1 


-173 


12.1 


14.4 


19 


-1.37 


PCII-6 


38.2 


25.5 


150.0 


-32.0 


2379 


2399 


60 


273 


2.37 


11.6 


2.7 


-93 


20.9 


15.4 


14 


-0.93 


PCII-7 


38.2 


25.5 


150.0 


-32.0 


2379 


2399 


60 


273 


2.37 


11.0 


2.3 


-66 


24.7 


14.5 


14 


-1.10 


PCil-8 


39.7 


28.2 


150.0 


-29.0 


2045 


2065 


59 


265 


2.30 


13.9 


6.8 


-57 


10.2 


12.8 


22 


-0.87 


PCII-9 


30.0 


0.0 


157.0 


-58.3 


2442 


2444 


173 


987 


4.16 


10.1 


2.4 


-177 


19.1 


15.9 


18 


-1.53 


PCII-IO'' 


163.8 


48.0 


162.4 


59.2 


2390 


2410 


150 


672 


3.86 


6.6 


2.5 


-132 


11.7 


4.6 


22 


-2.06 


PCII-11 


129.6 


53.9 


164.3 


37.2 


2316 


2331 


93 


425 


3.33 


10.0 


2.4 


-13 


10.9 


9.8 


20 


-1.04 


PCII-12 


124.5 


38.0 


183.4 


32.6 


2670 


2674 


83 


514 


3.46 


10.1 


1.5 


29 


19.0 


10.6 


16 


-0.90 


PCII-13 


64.8 


6.6 


187.0 


-29.5 


2805 


2826 


65 


353 


1.99 


11.9 


1.9 


20 


28.3 


16.1 


12 


-0.81 


PCII-14 


64.8 


6.6 


187.0 


-29.5 


2805 


2826 


65 


353 


1.99 


12.3 


1.7 


44 


31.3 


17.6 


14 


-0.82 


PCII-15 


116.9 


28.0 


192.4 


23.9 


2055 


2075 


35 


254 


3.12 


14.3 


1.1 


44 


9.1 


12.1 


10 


-0.73 


PCII-16 


139.4 


30.4 


195.6 


43.5 


2381 


2401 


114 


527 


3.50 


7.6 


2.0 


-103 


13.0 


10.1 


11 


-1.36 


PCII-17 


127.7 


24.4 


199.8 


32.0 


2315 


2330 


83 


431 


3.05 


12.5 


2.5 


-40 


18.1 


17.7 


10 


-1.46 


PCII-18 


165.6 


28.6 


203.1 


65.9 


2855 


2870 


151 


764 


3.21 


8.6 


4.7 


-157 


13.3 


14.2 


10 


-1.37 


PCII-19 


139.9 


22.2 


206.6 


41.9 


2304 


2319 


102 


609 


2.72 


14.5 


3.4 


-55 


27.3 


16.3 


11 


-1.16 


PCII-20 


132.6 


6.1 


221.5 


29.2 


2317 


2332 


69 


470 


3.04 


10.9 


1.6 


71 


13.2 


11.4 


15 


-0.96 


PCII-21'' 


134.0 


3.2 


225.2 


29.0 


2888 


2913 


74 


514 


3.50 


10.3 


1.4 


85 


14.9 


4.0 


17 


-1.02 



Note. — Kinematic and metallicity data for all class II peak detection ECHOS from [5o9 . The columns are: right ascension, declination, galacti( 
SEGUE faint plate number, number of SEGUE MPMSTO spectra along the indicated line of sight, number of photometric MPMSTO candidates along tl 
line of sight, median distance of the stars in the ECHOS, the error in that estimate, the radial velocity of the ECHOS, the velocity dispersion of the 
ECHOS, the number of MPMSTO stars kincmatically associated with the ECHOS, the iron metallicity of the ECHOS, the error in that estimate, th( 
the halo along the indicated line of sight, the error in that estimate, the Q-enhancement of the ECHOS, the error in that estimate, the ci-enhancement c 
indicated line of sight, and the error in that estimate. 

Monoccros Stream 

^ fGrillmair k. DionatosI 1(2003) Stream 



